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Channel  Diversity  in  Random  Wireless  Networks 

Kostas  Stamatiou,  Member,  IEEE,  John  G.  Proakis,  Life  Eellow,  IEEE,  and  James  R.  Zeidler,  Eellow,  IEEE 


Abstract — The  goal  of  this  paper  is  to  explore  the  benefits  of 
channel  diversity  in  wireless  ad  hoc  networks.  Our  model  is  that 
of  a  Poisson  point  process  of  transmitters,  each  with  a  receiver 
at  a  given  distance.  A  packet  is  divided  in  blocks  which  are 
transmitted  over  different  subbands  determined  by  random  fre¬ 
quency  hopping.  At  the  receiver,  a  maximum-likelihood  decoder 
is  employed  to  estimate  the  transmitted  packet/codeword.  We 
show  that,  if  L  is  the  Hamming  distance  of  the  error  correction 
code  and  e  is  a  constraint  on  the  packet  error  probability,  the 
transmission  capacity  of  the  network  is  proportional  to 
when  e  ^  0.  The  proportionality  constant  depends  on  the 
code  selection,  the  packet  length,  the  geometry  of  the  symbol 
constellation  and  the  number  of  receive  antennas.  This  result 
implies  that,  at  the  cost  of  a  moderate  decoding  complexity,  large 
gains  can  be  achieved  by  a  simple  Interference  randomization 
scheme  during  packet  transmission.  We  also  address  practical 
issues  such  as  channel  estimation  and  power  control.  We  find 
that  reliable  channel  information  can  be  obtained  at  the  receiver 
without  significant  rate  loss  and  demonstrate  that  channel- 
inversion  power  control  can  increase  the  transmission  capacity. 

Index  Terms — Frequency  hopping,  interference  diversity,  bit- 
interleaved  coded  modnlation  (BICM),  Poisson  point  process 


I.  Introduction 

THe  study  of  random  wireless  networks  has  recently 
gathered  a  lot  of  attention  in  the  research  community  [1]. 
The  main  theme  of  this  work  is  the  use  of  tools  from 
stochastic  geometry  in  order  to  characterize  the  performance 
of  an  ensemble  of  networks,  for  different  physical,  medium 
access  control  and  network  layer  strategies.  A  central  modeling 
assumption  is  that  the  network  consists  of  a  Poisson  point 
process  of  transmitters,  and  each  transmitter  (TX)  has  a 
corresponding  receiver  (RX)  at  a  given  distance.  A  popular 
metric  that  quantifies  the  network  performance  is  the  trans¬ 
mission  capacity,  defined  as  the  maximum  spatial  density  of 
transmissions,  multiplied  by  their  rate,  such  that  a  constraint 
on  the  packet  error  rate  is  satisfied  [2]. 

In  the  majority  of  existing  papers  (see  overview  in  [1]), 
interference  from  concurrent  transmissions  is  considered  as 
noise  and  an  outage  probability  approach  is  taken  to  model 
packet  successes:  Given  the  TX  locations  and  the  channels 
between  the  TXs  and  the  reference  RX,  which  are  assumed 
to  be  constant  during  the  transmission  of  a  packet,  a  packet  is 
successfully  received  if  the  signal-to-interference-ratio  (SIR) 
is  larger  than  a  certain  threshold.  In  information-theoretical 
terms,  assuming  the  TXs  are  sending  symbols  from  a  Gaussian 
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alphabet,  a  packet  reception  occurs  when  the  channel  mutual 
information  is  at  least  equal  to  the  desired  information  rate  [3]. 

In  this  paper,  as  in  [4],  we  take  an  alternative  approach 
and  explore  the  impact  of  channel  randomization  within  the 
transmission  of  a  packet.  Our  motivation  stems  from  the  well 
known  fact  that  channel  diversity  can  be  exploited  through 
error  correction  coding  in  order  to  yield  performance  gains. 
Specifically,  we  consider  a  physical-layer  scheme  based  on 
random  frequency  hopping  (FH),  bit-interleaved  coded  mod¬ 
ulation  (BICM)  [5],  and  maximal  ratio  combining  (MRC)  to 
exploit  spatial  diversity  at  the  RX.  FH  is  preferred  over  Direct 
Sequence  (DS)  Spread  Spectrum  (SS)  as  the  multiple-access 
(MA)  scheme  primarily  because  DS-SS  suffers  from  the  near- 
far  problem  in  the  decentralized  environment  of  an  ad  hoc 
network  [6].  Moreover,  it  is  known  that  FH  combined  with 
coding  exploits  frequency  diversity,  if  the  hopping  distance  is 
larger  than  the  coherence  bandwidth  of  the  channel  fading, 
and  interference  diversity,  as  the  set  of  interfering  TXs  over 
each  dwell  is  potentially  different.  FH  also  simplifies  the  MA 
protocol,  as  no  coordination  is  required  between  different  TX- 
RX  links.  Regarding  the  choice  of  coding  and  modulation 
strategy,  BICM  is  spectrally  efficient  and  well  suited  for  fading 
and  interference  channels,  as  the  diversity  order  is  determined 
by  the  Hamming  distance  of  the  employed  convolutional 
code  [5].  Finally,  MRC  only  requires  channel  state  information 
at  the  RX,  which  is  desirable  given  that  the  transmission 
subband  changes  frequently  across  time. 

We  analyze  the  performance  of  the  aforementioned  scheme 
in  terms  of  the  codeword/packet  error  probability  and  evaluate 
the  transmission  capacity.  Since  an  averaging  over  different 
channel  states  takes  place  within  a  packet,  the  information  rate 
of  the  typical  TX-RX  link  is  upper-bounded  by  the  ergodic 
capacity  for  which  we  provide  tight  upper  and  lower  bounds. 
We  also  address  practical  physical-layer  issues  such  as  channel 
estimation,  power  control  (PC)  and  channel  correlation,  and 
assess  their  effect  on  the  performance  via  simulation. 

A.  Related  work 

Several  papers  have  dealt  with  the  performance  analysis  of 
coded  FH  systems  under  MA  interference,  as  well  as  partial- 
band  interference  (see,  e.g.,  [7],  [8]  and  Chapter  12  of  [9]). 
A  common  feature  of  such  systems  is  that  the  performance 
can  be  dramatically  improved  if  the  decoder  is  aware  of  the 
interference  levels  across  the  codeword.  If  a  Reed-Solomon 
(RS)  code  is  employed,  the  decoder  declares  an  erasure  when 
a  symbol  has  been  “hit”;  in  the  case  of  soft-decoding,  the 
metrics  in  the  Viterbi  decoder  are  weighted  by  the  respective 
SIRs.  In  [10],  RS  coding  combined  with  FSK  modulation 
is  considered  in  a  Poisson  field  of  interferers  and  the  im¬ 
pact  of  the  code  rate  on  the  information  efficiency,  i.e.,  the 
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product  (packet  success  probability)  x  (transmission  distance) 
X  (rate),  is  explored.  More  recently,  [11]  has  extended  the 
work  in  [10]  to  accommodate  differential  unitary  space-time 
modulation  and  unknown  fast  time-varying  channels. 

The  use  of  spread-spectrum  (SS)  communication  for  ad  hoc 
networks  is  discussed  in  [6].  The  authors  make  an  argument 
against  interference  averaging  which  they  define  as  “. . .  using 
DS-SS  or  fast  FH^  to  proportionally  reduce  the  interference 
level”  and  advocate  hopping  at  the  packet  level,  or  interference 
avoidance  (lA),  as  the  preferable  MA  scheme  for  ad  hoc  net¬ 
works.  While  the  near-far  problem  of  DS-SS  in  a  decentralized 
environment  is  clear,  it  is  not  obvious  why  slow  FH  might 
be  preferable  to  fast  FH,  apart  from  the  fact  that  slow  FH 
potentially  induces  less  overhead  in  terms  of  code  acquisition 
and  synchronization. 

B.  Contributions 

This  paper  demonstrates  that  considerable  gains  in  terms 
of  network  capacity  are  possible  by  combining  FH  during 
packet  transmission  and  error  correction  coding  of  modest 
complexity.  If  L  is  the  Hamming  distance  of  the  convolutional 
code  employed  at  the  TX,  A  is  the  density  of  TXs  and  M  is 
the  number  of  subbands,  we  show  that,  as  A/M  ^  0,  the 
codeword  error  probability  follows  the  power  law 
T]  >  0.  This  leads  to  the  conclusion  that,  for  e  ^  0,  where 
e  is  the  constraint  placed  on  the  codeword  error  probability, 
the  transmission  capacity  is  proportional  to  The  propor¬ 
tionality  constant  depends  on  the  selected  code,  the  codeword 
length,  the  geometry  of  the  symbol  constellation,  as  well  as 
the  term  where  N  is  the  number  of  RX  antennas  and 

6  >  2  is  the  propagation  exponent.  We  also  derive  upper  and 
lower  bounds  on  the  ergodic  capacity  C  of  the  typical  TX- 
RX  link.  Specifically,  we  show  that  C  >  |  log2  (/iA^^/^^), 
where  /r  >  0  is  an  appropriately  defined  constant. 

Practical  physical-layer  issues  are  discussed  such  as  chan¬ 
nel  estimation,  PC  and  channel  correlation.  We  demonstrate 
via  simulation  that,  with  an  acceptable  rate  loss  due  to  the 
transmission  of  pilot  symbols,  the  channel  fading  and  the 
interference  power  can  be  accurately  estimated  for  decoding. 
With  respect  to  PC,  it  is  shown  that  channel  inversion  can 
actually  enhance  the  capacity,  since  the  error  correction  code 
protects  the  RX  from  the  deep  fades  of  its  nearby  interferers. 
Finally,  the  impact  of  the  channel  correlation  is  assessed  as  the 
number  of  subbands  and/or  the  number  of  dwells  is  decreased 
and  it  is  shown  that  the  gains  compared  to  slow  FH  are  still 
significant. 

C.  Paper  organization  and  notation 

The  rest  of  the  paper  is  organized  as  follows.  Section  II 
introduces  our  system  model  in  detail.  In  Section  III  we  derive 
the  statistics  of  the  SIR  and  determine  the  performance  of 
the  decoder  under  perfect  CSI.  The  transmission  capacity  is 
defined  and  evaluated  in  Section  IV.  Section  V  discusses 
practical  physical-layer  considerations  and  Section  VI  presents 
our  numerical  results.  Section  VII  concludes  the  paper. 

’Fast  FH  refers  to  hopping  on  the  order  of  a  symbol  or  a  few  symbols, 
while  slow  FH,  or  interference  avoidance,  refers  to  hopping  at  the  packet 
level. 


TABLE  I 

Commonly  used  symbols 


Symbol 

Meaning 

R 

Distance  of  typical  TX-RX  link 

A 

Density  of  TXs 

M 

Number  of  subbands 

N 

Number  of  RX  antennas 

b 

Path-loss  exponent  (b  >  2) 

a  =  2/b 

Stability  exponent 

Lb 

Number  of  information  bits  per  packet 

Me 

Number  of  bits  per  constellation  symbol 

Rc 

Rate  of  convolutional  code 

Td 

Number  of  data  symbols  per  dwell 

Tv 

Number  of  pilot  symbols  per  dwell 

D 

Number  of  dwells 

L 

Hamming  distance  of  convolutional  code 

Pi 

Probability  of  length-/  error  event 

Wl 

Number  of  length-/  error  events 

Pe 

Probability  of  packet  error 

C 

Ergodic  capacity  of  typical  TX-RX  link 

Te 

Transmission  capacity  under  constraint  Pe  =  e 

A  real  (circularly  symmetric  complex)  Gaussian  random 
variable  (r.v.)  x  with  mean  0  and  variance  is  denoted  as 
X  ~  A/'(0,  cr^)  {x  ~  CA/'(0,  (T^)).  A  central  chi-square  r.v.  x 
with  parameter  1/2  and  n  degrees  of  freedom  is  denoted  as 
X  ~  I„  is  the  n  X  n  identity  matrix.  (•)"'"  and  (•)^  denote 
the  transpose  and  conjugate  transpose  operations,  respectively. 
[X]„  t  denotes  the  (n,  t)  element  of  matrix  X  and  \x\n  denotes 
the  element  of  vector  x.  The  symbol  is  employed  to 
denote  asymptotic  equality  of  two  functions.  A  list  of  symbols 
commonly  used  throughout  the  paper  is  provided  in  Table  I. 

II.  System  model 

A.  General 

We  consider  a  network  of  TXs,  each  with  a  RX  at  a  fixed 
distance  R  and  random  orientation.  The  locations  of  the  TXs 
are  determined  according  to  a  homogeneous  Poisson  point 
process  H  =  {zi}  C  z  €  N,  of  density  A.  The  TXs  send 
packets  to  their  corresponding  receivers  concurrently  and  in  a 
synchronized  manner.  Typically,  the  locations  of  the  nodes  are 
constant  for  at  least  the  duration  of  a  packet. 

The  bandwidth  is  divided  into  M  subbands.  The  channel 
between  a  typical  TX-RX  pair  over  a  subband  comprises  flat 
Rayleigh  fading  and  path  loss  according  to  the  law  r~^,  where 
6  >  2  is  the  propagation  exponent^.  We  assume  that  the 
coherence  bandwidth  of  the  fading  is  equal  to  the  width  of 
a  subband,  while  the  coherence  time  is  at  least  equal  to  the 
duration  of  a  dwell,  which  will  be  defined  shortly.  We  also 
consider  an  interference-limited  scenario,  i.e.,  the  TX  power  is 
sufficiently  large  that  additive  noise  at  the  RX  can  be  ignored. 
Initially,  we  assume  that  the  TX  power  is  the  same  for  all 
TXs.  Issues  of  PC  to  compensate  for  long-term  fading,  e.g., 
shadowing,  are  discussed  in  Section  V. 

Assume  that  a  packet  corresponds  to  Tb  binary  information 
bits,  which  are  the  input  to  a  convolutional 

^The  constraint  b  >  2  is  required  in  order  for  the  interference  power  to  be 
finite  [1],  [12], 
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encoder  of  rate  Rc  <  1^.  The  bits  ci, . . .  ,Ck'  ,■  ■  ■ 
of  the  output  codeword  are  interleaved  and  Gray-mapped  to 
symbols  xi, . . . ,  Xk,  •  ■  • ,  XL^/(^fi^Mc)  from  a  complex  PSK  or 
QAM  constellation  X  of  size  \X\  =  2^“=,  zero  mean  and 
unit  average  power.  We  assume  that  the  one-to-one  interleaver 
mapping  k'  ^  {k,  jk'),  where  jk'  =  1,  ■  •  ■ ,  Me  is  the  position 
of  Ck'  in  the  symbol  Xk,  is  known  at  the  RX.  Next  in  the  TX 
chain,  the  symbol  sequence  is  divided  in  D  =  L^/ {RcMcTd) 
groups"^  of  size  Td  and  each  group  is  transmitted  in  a  dwell, 
over  a  subband  which  is  randomly^  selected  with  probability 
1/M.  If  we  denote  the  data  symbols  of  the  dwell, 
d  =  as  xj,  the  sequence  constitutes 

a  packet.  For  convenience,  we  assume  all  transmissions  are 
synchronized  at  the  dwell  level  (the  issue  of  asynchronous 
transmissions  is  discussed  in  Section  V). 

Consider  the  typical  TX,  located  at  zq,  and  its  corresponding 
RX,  both  specified  by  index  0,  i.e.,  TXq  and  RXq.  If  RXq  is 
equipped  with  an  antenna  array  of  size  N  >  1,  the  received 
data  matrix  in  dwell  d  is® 

Yd,o  =  hd, oxj^o  +  X!  (^d,ir~^^'^e''^*hd,ixj^,  (1) 

i#0 

where  hd.o  ~  CAf(0,lAr)  is  the  fading  vector  between  TXq 
and  RXo;  Cd.i  is  the  indicator  of  the  event  that  the  TX  located 
at  Zi  (denoted  as  TXQ  and  at  distance  from  RXq,  occupies 
the  same  subband  as  TXq  in  dwell  d,  i.e.,  iP{ed,i  =  1)  =  1/M 
and  P(ed,i  =  0)  =  1  —  1/M;  hd.i  ~  CN{Q,1n)  is  the 
fading  vector  between  TX^  and  RXq;  £cJ ^  is  the  group  of  data 
symbols  transmitted  by  TX^  in  dwell  d;  and  c/i  is  a  random 
phase,  uniformly  distributed  in  [0, 27r),  which  models  the 
phase  offset  between  the  RXq  and  TX^.  Note  that  the  subscript 
d  in  the  fading  vectors  indicates  that,  in  general,  depending  on 
the  coherence  time,  these  may  vary  independently  from  dwell 
to  dwell. 

Let  Wrf  o  denote  the  interference  term  in  (1).  Since  the 
elements  of  Xd.i  are  zero-mean  and  independent,  the  same 
holds  for  the  elements  of  Moreover,  provided  that 

[xdd]t  is  selected  from  a  PSK  constellation  of  unit  power, 
[hd,/\n[xdAt  ~  CA/'(0, 1),  n  =  1, . . . ,  N ,  t  =  1, . . .  ,Td. 
Given  {ed,i,ri},  it  follows  that  [Wd,o]n.i  ~  CAf{0,Id,o), 
where  Id.o  —  ^d,iT//^  is  the  interference  power  seen 

by  RXq  in  dwell  d.  In  the  case  of  QAM  constellations,  we 
approximate  [W^  o]n,t  as  Gaussian  for  the  sake  of  tractability. 

RXq  can  obtain  knowledge  of  hd.o  and  Id.o  with  the  help 
of  pilot  symbols  which  are  transmitted  at  the  beginning  of  the 
dwell.  Presently,  we  assume  that  they  are  perfectly  known;  a 
straightforward  channel  estimation  algorithm  is  presented  in 
Section  V. 


B.  Equivalent  channel  model  and  decoding 

The  reference  RX  performs  MRC,  i.e.,  it  evaluates  the 
product  Yd,o-  From  (1),  we  have 


\hd 


^Yd  =  Xd 


K, 


\hd 


(2) 


where  we  have  omitted  the  index  0  in  order  to  simplify  the 
notation.  Denoting  the  f**'  columns  of  and  as  y^.  and 
Wt,  respectively,  we  have  the  following  equivalent  channel 
model  for  data  symbol  Xk,  k  =  {d—  l)Td  -f  1, . . . ,  dTd,  which 
is  transmitted  in  dwell  d 


Vk  =  Xk  +  Wk,  (3) 

where  yk  =  ~ 

CA/’(0, /rflAf),  it  follows  that  Wk  ~  cjV(0,7/"^),  given  the 
equivalent  SIR  7^  =  ,  where,  by  definition,  d  =  [^1- 

The  r.v.  ad  =  ||^d|P  is  chi-square  distributed  with  2N  degrees 
of  freedom,  i.e.,  ad  ~  X2N-  addition,  due  to  the  fact  that 
the  locations  of  the  inteferers  in  each  dwell  are  a  realization 
of  a  Poisson  point  process  with  density  A/M,  it  is  known  that 
Id  is  an  a-stable  r.v.  with  stability  exponent  a  =  2/6  [12], 
[14].  Its  moment  generating  function  (mgf)  is  [12] 

$/(s)  =  ]E[e“'^®]  =  ,  s  >  0,  (4) 


where 

6  4  7rr(l  -  a)i?2  (5) 


and  r((^),  C  >  denotes  the  gamma  function.  The  sequence 
{yk,lk}k=i  input  to  the  decoder  which  decides  that 

the  codeword  c  was  transmitted  according  to  the  simplified 
maximum-likelihood  (ML)  criterion  [5,  Eq.  (9)] 


c  =  arg  min 


ib/fle 

E 

fc'=i 


7fc  min  {|yfe-a:p},  (6) 


x&X, 


where  X//l  denotes  the  set  of  constellation  symbols  that  have 
bit  Ck'  at  position  jk',  jk'  =  1,  •  •  • ,  Me.  The  weighting  of  each 
distance  metric  by  the  respective  SIR  reflects  the  confidence 
of  the  decoder  in  that  metric.  Note  that  the  summation  is  over 
the  bit  index  k' ,  which  corresponds  to  a  unique  pair  {k,jk') 
through  the  interleaver  mapping. 


III.  Analysis 


^We  assume  that  the  encoder  is  trellis-terminated  [13],  i.e.,  it  is  forced  to 
Start  from  and  end  at  the  zero-state.  This  results  in  a  small  rate  loss  which  is 
not  taken  into  account. 

“^For  convenience,  we  assume  that  L\^/ (RcMc)  and  D  are  integer. 

^In  reality,  the  hopping  pattern  is  determined  pseudorandomly  and  is  known 
at  the  RX.  The  model  of  random  FH  is  convenient  for  analytical  purposes. 

^We  assume  that  the  average  received  power  R~^  per  antenna  is  known  at 
RXo .  We  have  taken  it  into  account  in  the  interference  portion  of  the  received 
signal  because  it  is  convenient  in  terms  of  notation. 


This  section  is  devoted  to  the  performance  analysis  of  de¬ 
coder  (6).  In  order  to  round  out  the  analysis,  in  Section  III-D, 
we  also  derive  upper  and  lower  bounds  to  the  ergodic  capacity 
of  channel  (3).  Since  the  scheme  presented  in  Section  II 
induces  an  “averaging”  over  different  channel  states  within 
the  packet,  the  ergodic  capacity  is  an  upper  bound  to  the 
information  rate  of  the  typical  TX-RX  link. 
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A.  Decoder  performance 

The  codeword  (or  frame)  error  probability  (FEP)  of  decoder 
(6),  Pe,  is  upper-bounded  as  [13] 

Pe<U  Y. 

l=LX+l,... 


where  Pi  is  the  probability  of  a  length-^  error  event,  or  pairwise 
error  probability,  and  wi  is  the  number  of  length-!  error  events. 
The  minimum  length  of  an  error  event  L,  i.e.,  the  Hamming 
distance,  as  well  as  the  weight  distribution  {w;}  depend  on 
the  particular  code  employed. 

We  now  assume  that,  due  to  random  interleaving,  the 
sequence  of  I  symbols  that  correspond  to  a  sequence  of  I  coded 
bits  encounter  independent  SIR  conditions  (this  assumption  is 
discussed  in  Section  V).  From  [5],  [15],  Pi  can  be  upper- 
bounded  by 


TT 


E 

{x,x')G‘^ 


\  4  sin^  ^  j 


d9,  (8) 


where  <i>.y(s)  =  E[e“'*'®],  s  >  0.  (The  time  index  k  has 
been  removed  as  it  is  of  no  consequence.)  {x,x')  are  all 
possible  Mc2^‘=  nearest-neighbor  pairs  in  X  which  have 
complementary  bits  in  position  j,  j  =  1, . . . ,  Me,  e.g.,  for 
QPSK  constellations  with  Gray  mapping,  all  such  pairs  are 
at  the  minimum  constellation  distance  dmin-  Note  that,  in  the 
limit  of  large  £^[7],  the  bound  is  tight,  i.e..  Pi  ~  Pi  [5]. 

From  (7),  a  further  upper  bound  to  Pe  is 

Fe  =  £b  E 

l=L,L+l,... 

The  evaluation  of  Pi  in  (8)  is  possible  by  numerical  integra¬ 
tion,  provided  that  is  known.  The  derivation  of 

is  the  topic  of  the  following  subsection. 


B.  Statistics  of 

The  probability  density  function  (pdf)  of  7  is  given  by  the 
following  lemma. 


Lemma  1  The  pdf  of  the  SIR  7  is 


f  \P^\  (XS 

^  '  n—1  ^ 


(10) 


a;  >  0,  where 

n 

=E(-l)^ 


m—1 
A 


n 

m 


{am)N,  n=l,...,N  (11) 


and  {am)]\[  =  am  . . .  {am  —  N  +  1)  is  the  falling  sequential 
product. 

Proof:  The  cumulative  distribution  function  (cdf)  of  7, 
F-f{x),  is  by  definition 

^4-00 

Fj{x)  =  P(a  <xl)=  /  Fa{xy)fi{y)dy,  (12) 
Jo 


where 

Fa{x)  =  1  -  e 


V  — =  i-EEE 

nl  {N  —  1)!  ’ 


x>0  (13) 


n—0 


is  the  cdf  of  a  and  r(C,  a;)  =  e  ^df,  a;  >  0,  i 


IS 


the  incomplete  gamma  function.  Substituting  (13)  in  (12)  and 
taking  the  derivative,  the  pdf  of  7  is  given  by 

1  f+^dr{N,xy) 


fiix)  =  - 


{N-l)Uo 


da; 


-//(y)dy 


^4-00 


(N-iy.Jo 

{—l)^x^~^  d^^i{x) 
{N-l)l  da;^ 


fi{y)y  e  ®^d?/ 


where  we  have  used  the  identity 
p.  951]  and  the  Laplace  transform  property' 


dr(C,a^)  ^ 


dx 


(14) 


[16, 


fi{y)y 


^  ^  ds^  ■ 


From  the  identity  for  the  derivative  of  a  composite 

function  [16,  Eq.  (0.430.1),  p.  24],  after  some  algebra  we 
obtain 


^  n\  \M 

n—1  ^ 


ds^ 


with  (3^  given  by  (11).  From  (15)  and  (14) 


Md  = 


1 


e 


N 

■E' 

n—1 


(15) 


(16) 


(W-1)!  a: 

In  order  to  derive  (10)  from  (16),  we  need  to  show  that 
{—1)^P^  >  0.  Once  again,  using  the  identity  for  the 
derivative  of  a  composite  function,  can  be  written  as  the 
following  derivative  evaluated  at  a;  =  1: 

d'^  (1  —  x°'y 


= 


dx^ 


(17) 


From  (17),  the  following  iterative  relation  can  be  proved  for 
N>2 


N 


=  E 


mi  =  l 


N 

mi 


mi  oN—mi 


n—1 


(18) 


By  successive  application  of  (18),  we  obtain 

(-1)^/3^ 


Nl 


N 

E 


N —mn  —  2 


E 


where  to„  =  TV  — m„_i - -  —  mi.  Due  to  (11),  (— = 

(— 1)^4-1q,(q,  _  1)  _  (q,  _  _|_  1)^  Since  a  =  2/6  <  1,  it 

follows  that  {—l)^py  >  0.  Therefore,  (—1)'^/?^  >  0  for 
n  =  l,...,TV.  ■ 

As  expected,  increasing  the  spatial  diversity  order  N  in¬ 
creases  the  value  of  f-y{x),  as  more  positive  terms  are  added 
to  the  polynomial  in  (10).  Eq.  (10)  enables  the  numerical 
evaluation  of  4)-y(s)  using  Gauss-Laguerre  quadrature  for  any 
s  >  0. 


^This  identity  is  also  employed  in  [17],  in  order  to  derive  the  cedf  of  7. 
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C.  Approximations 

The  numerical  evaluation  of  Pi  from  (8)  does  not  provide 
insight  on  how  the  decoder  performance  depends  on  the 
system  parameters.  In  this  section,  we  examine  the  decoder 
performance  when  A/M  — >  0,  i.e.,  the  interferer  point  process 
in  each  dwell  is  sparse.  This  implies  that  the  network  is 
operated  in  a  regime  of  small  FEP,  i.e.,  typically  //=  <  e, 
with  e  <  0.1. 

Let  i3(Ci,C2),  Ci)C2  >  0  denote  the  beta  function.  Our 
main  result  is  stated  in  the  following  proposition. 


Proposition  1  If  X/M  ^  0,  then  Pi  = 
where  Pi  is  defined  in  (8)  and  rj  is  a  positive  constant. 
Moreover, 


■{( 


M  ) 


/  I  I 

Pi  — - B  [  od  -\-  —,al  — 


where 


xB{N  —  a,  a)  X-kE? 


dl 


M 


(19) 


4  = 


E 


Mc2^«  \x  —  X' 

{x,x')^X 


J\2ol 


(20) 


Proof:  Omitting  the  term  e  in  (10),  an  upper  bound 

to  <h7(s)  for  all  s  >  0  is 


^  1 

n\  \M 

''  '  n—1 


^an-l^-xs 


dx 


E 


14 


T(N)  ^  n! 

^  '  n—1 


aB{N  —  a,a)  /  X6 

r(i-a)  [m 
1  V  141 

r(iV)^2 


T(an). 


(21) 


Eq.  (21)  is  derived  by  successively  applying  the  identity  r(C  + 
1)  =  Cr(C)  to  show  that 


|/3f  I  =a{l-a)...{N-l-a)  = 


r(l-a) 

and  the  definition  of  the  beta  function 

T{N-a)T{a) 


r{N-a) 


B{N  —  a,a)  = 


m) 


(22) 


where  T{N)  =  {N  —  1)!.  Erom  (21)  and  (8),  we  can  see  that 
+  ®  for  A/M  — >  0,  with  p  appropriately 

defined. 

Note  that  the  bound  in  (21)  is  tight  as  A/M  ^  0.  Ignoring 
the  higher  order  terms  and  recalling  the  definition  of  <5  in  (5) 

A 


<I).y(s)  ~  ‘h7(s)  —  TraR^B{N  —  a,  a)s 


M' 


(23) 


Substituting  (23)  in  (8) 


Pi  ~  (  TraR^B{N  -  a,  )  — 


\  \  I  ,al  /•ir/2 

-2  '  I  (sin  61)2“'  d6», 


M 


TT 


where  is  defined  in  (20).  Employing  the  identity  [16, 
p.  412] 


ir/2 


(sin6»)2“'  d6» 


22“'-1b(M  +  1/2,M  +  1/2) 


we  obtain  (19).  ■ 

Remarks:  Proposition  1  states  that  for  A/M  — >  0  the 
Hamming  distance  of  the  convolutional  code  determines  the 
diversity  order,  i.e.,  the  slope  of  the  curve  of  P^  vs.  A/M  in 
log-log  coordinates.  This  result  is  reminiscent  of  the  perfor¬ 
mance  of  BICM  in  a  Rayleigh  fading  channel  with  AWGN. 
The  difference  is  that,  in  the  context  of  an  inteference-limited 
random  network,  the  average  signal-to-noise  ratio  is  replaced 
by  where  iVeff  =  44  expected  number  of  TXs 

in  the  transmission  range  per  subband. 

Eq.  (19)  also  reveals  that  MRC  introduces  an  array  or 
coding  gain  through  the  factor  B{N  —  a,a)^.  To  obtain 
further  insight,  we  examine  the  trend  of  the  beta  function  for 
large  N.  By  Stirling’s  approximation,  for  large  ^  we  have 
r(C)  ~  therefore 


T{N-a) 

m) 


However,  it  is  easy  to  verify  that  limAr^oo  (1  ~  ^  = 

e““,  so,  from  (22),  B{N  —  a,a)  ~  r(Q;)A^““.  As  a  result, 
for  large  N,  the  coding  gain  is  proportional  to 

A  final  observation  is  that,  similarly  to  [5],  the  parameter 
is  the  harmonic  mean  of  the  minimum  squared  Euclidean 
distance  between  the  nearest-neighbor  pairs  defined  in  Sec¬ 
tion  III- A,  raised  to  the  stability  exponent  a. 


D.  Ergodic  capacity 

The  ergodic  capacity  of  (3)  is 

^-1-00 

c  =  /  /7(a^)log2(l  +  a;)da;.  (24) 

Jo 

A  closed  form  expression  appears  hard  to  obtain  due  to  the 
complicated  nature  of  (10);  nevertheless,  the  integral  can  be 
evaluated  numerically  with  Gauss-Laguerre  quadrature.  The 
following  proposition  provides  upper  and  lower  bounds  on  C. 


Proposition  2  The  ergodic  capacity  of  (3)  is  upper-bounded 
as 


C<c  =  log2 


l^i-f  A^r 


(25) 


and  lower  bounded  as 


C>C=-  log, 
a 


Hn-1 

^  In  2 


(26) 


where  T  =  0.577  ...  is  the  Euler-Mascheroni  constant  and 


Hr,  = 


^  k 

k=l 


n>  1 


(27) 


n  =  0 
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is  the  harmonic  number.  This  bound  is  tight,  i.e.,  C  ~  C_ 
for  A  — >  0.  A  looser  lower  bound  is 


C  = 


7V“MA 
e^A(5  / 


(28) 


Proof:  The  upper  bound  is  derived  by  noting  that  ]E[7]  = 
E[a]E  [i]  =  NE  [i],  so,  for  TV  =  1,  E  [i]  =  £[7],  Setting 
=  1  in  (10),  we  have 


E 


1 

7 


X6a 

~w 


r+00 


da; 


By  applying  Jensen’s  inequality*  on  (24),  we  obtain  (25). 

For  the  derivation  of  the  lower  bound,  we  employ  the 
inequality  ln(l  +  a;)  >  In  a;  for  x  >  0.  Combining  (14)  and 
(24)  and  integrating  by  parts,  we  obtain 


p  +  ca 

((TV- l)!ln2)C  =  (-1)^  /  a;^-i  ($/(a;))^^hna;dx 

Jo 

=  (-1)^  H  {x^-^  Ina;) 

k^O 

—  [  (4)7(3;))^^^  (a;'^“^  Ina;)^^  dx. 

Jo 


(fc)]  +°° 

0 


(29) 


For  convenience,  the  notation  (^/(x))^"^  has  been  adopted 
to  denote  the  if^  derivative  of  4>7(x)  and  likewise  for  the 
other  functions.  After  some  algebra,  we  can  show  that,  for 
A:  =  1,...,TV-  1,  >  1, 


(x^-i  Inx)^'"^  =  {N-1)...{N-  Inx 

^x^v-fe-i  ^  (TV-Ti)...(TV-/fc_i),  (30) 


with  the  summation  taken  over  all  permutations  of  the  vector 
(Ti, . . . ,  Ik-i),  Ij  =  j  =  -  1.  When  k  = 

A^-  1 

(x^-^  lnx)^^“^^  ={N-  l)!(lnx  +  Hn_i),  (31) 


where  iT„  is  defined  in  (27).  From  (15)  and  (30),  we  can  show 
that  the  first  term  in  (29)  is  zero.  Hence 

7+00 


C  = 


X6a 

M\n2 


x“  ““  Inx  dx  + 


Hn-1 

In  2 


from  which  (26)  follows  by  use  of  [16,  Eq.  (4.331.1),  p.  602]. 
Since  the  harmonic  number  is  lower  bounded  as  [19] 

i7jv>lnTV  +  r  + 

for  TV  >  2,  we  have  that 

Hn-1  =  Hn  ~  >lnTV  +  r  + 


2(Af  +  l) 
1 


2(W+1)  N 


The  latter  inequality  holds  because,  for  TV  >  2,  ^  —  2{n+i)  ^ 
2  <  F,  VTV  >  2.  As  a  result,  a  looser  lower  bound  to  the 
capacity  is  (28).  ■ 

Eq.  (28)  shows  that  C  is  a  linear  function  of  log2  (^),  with 
slope  \/ a  =  6/2  and  a  constant  term  log2  N  —  ^  log2  e. 


*  Jensen’s  inequality  was  also  employed  in  [18]  in  order  to  derive  an  upper 
bound  to  the  ergodic  capacity,  albeit  in  a  slightly  different  context. 


IV.  Network  metrics 

Having  evaluated  the  performance  at  the  link  level,  we 
now  turn  our  attention  to  network-wide  metrics.  The  network 
thoughput  per  unit  area  is  defined  as  the  spatial  density  of 
successful  transmissions  multiplied  by  their  rate  RcMc,  i.e., 

r  =  A(1  -  Pe)RcM,. 

Similarly  to  [2],  the  transmission  capacity  is  defined  as 
the  maximum  network  throughput  per  unit  area  such  that  a 
constraint  =  e  is  satisfied,  i.e.. 


T,  =  X,{1  -  €)RcM„  (32) 


where  A^  is  the  maximum  contention  density.  A  closed-form 
expression  for  A^  may  be  obtained  by  noting  that,  for  e  — >  0, 
the  Hamming  distance  error  events  dominate  the  EEP,  i.e., 
Pe  —  77  —  LbWLPL-  From  Proposition  1,  we  find  that 


/  e  \i/^  dj^M 

\k)  16“q;P(TV  —  Of,  a)7rP2  ’ 


(33) 


where 

K  =  ^J^B{aL  + 1/2, aL+ 1/2).  (34) 

27r 

The  maximum  contention  density  is  therefore  proportional  to 
as  e  ^  0.  This  result  is  a  manifestation  of  the  channel 
diversity  harnessed  through  FH  and  coding  at  the  expense  of 
spectrum.  As  seen  from  (33)  and  (34),  the  proportionality 
constant  depends,  among  other  system  parameters,  on  the 
selected  code,  i.e.,  L  and  wl,  the  codeword  length  Lb,  the 
geometry  of  the  symbol  constellation  and  the  number  of  RX 
antennas  TV.  Since  B{N  —  a,a)  «  TV““  for  increasing  TV, 
it  also  follows  that  Ae  is  proportional  to  TV“,  which  is  in 
agreement  with  the  scaling  law  derived  in  [17]. 

As  a  final  note,  the  network  throughput  is  upper-bounded 
by  f  =  AC,  where  C  is  the  link  ergodic  capacity  given  by 
(24).  Bounds  on  f  may  be  obtained  by  employing  the  results 
of  Proposition  2. 


V.  Physical-layer  considerations 

This  section  discusses  various  physical-layer  issues  with 
respect  to  the  system  model  presented  in  Section  II.  The 
influence  of  these  on  the  decoder  performance  is  assessed  via 
simulation  in  Section  VI. 


A.  Channel  estimation 

In  this  subsection,  we  discuss  the  important  issue  of  how 
the  decoder  obtains  estimates  of  hd  and  Id  in  the  dwell. 
Assume  that  a  header  of  Tp  pilot  symbols,  selected  from  a 
complex  PSK  constellation  of  zero  mean  and  unit  power,  and 
known  at  the  RX,  is  transmitted  at  the  beginning  of  the  dwell. 
If  this  header  is  denoted  as  p^,  then,  similarly  to  (I),  the 
received  pilot  matrix  is 

Yd  =  hdpZ  +  Wd 

where  Wd  =  ed,ir~''^‘^P^'hddPJd  [Wd]„,t, 

n  =  1,...,TV,  t  =  l,...,Tp  are  i.i.d.  with  [Wd]„,t  ~ 
CAf{0,Id),  given  Id. 
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An  estimate  of  Id  is 


h  = 


N{Tp  -  1) 


It  -  T^PdPd  )  Yd 

J-  r> 


obtained  by  finding  jointly  the  ML  estimators  of  hd  and  Id  [20, 
p.  182]  and  multiplying  the  latter  by  the  factor  Tp/{Tp  —  1) 
in  order  to  remove  the  bias.  Intuitively,  the  estimate  of  the 
interference  power  is  obtained  by  projecting  the  received 
signal  onto  the  subspace  which  is  orthogonal  to  the  pilot 
data  p.  Assuming  that  the  estimate  of  Id  is  accurate,  i.e.. 
Id  ~  Id,  the  minimum  mean-square  error  (MMSE)  estimate 
of  hd  is  [20,  p.  391] 


hd 


1 

Id  +  Tp 


^dp*d. 


Defining  Ed  =  hd  —  hd,  it  holds  that  hd  ~  CM  (|o, 

Ed  ~  CM  ^0,  and  hd,  Ed  are  independent.  We  can 

see  that  the  estimate  of  hd  is  accurate  if  Tp  ^  Id- 

Following  the  estimation  of  hd  and  Id,  the  RX  performs 

MRC,  i.e.,  it  evaluates  the  product  -  Y^.  From  (1),  we 

W^d  II 

have 


h^ 


h 


Lumping  the  channel  estimation  error  term  with  the  interfer¬ 
ence,  we  can  show  that  the  equivalent  channel  model  follows 
(3),  with  the  SIR  defined  as 


ad 


^k,csi  —  r 

dd 


C.  Synchronization 

The  assumption  of  synchronization  between  different  TX- 
RX  pairs  in  the  network  facilitates  the  analysis,  as  the  interfer¬ 
ence  power  is  the  same  throughout  each  dwell.  We  anticipate 
that  the  results  of  our  analysis  will  also  hold  in  the  case 
of  asynchronous  transmissions,  provided  that  the  interference 
power  can  be  estimated  accurately.  This  may  necessitate  the 
insertion  of  more  pilot  symbols  throughout  the  dwell  and, 
consequently,  the  loss  of  rate.  The  study  of  these  issues  is 
beyond  the  scope  of  the  present  paper. 

D.  Power  control 

In  [2],  it  was  shown  that  channel-inversion  PC  increases  the 
outage  probability  in  a  random  network  with  fading.  In  order 
to  examine  the  effect  of  PC  in  the  context  of  this  paper,  we 
assume  that  the  typical  TX-RX  link  is  subject  to  lognormal 
shadowing,  i.e.,  the  transmitted  signal  is  multiplied  by  a  r.v. 
S  =  where  ^  ~  ^(0, 1)  and,  typically,  Cs  =  6  — 

8  dB.  Note  that  S  (like  the  node  locations)  is  assumed  to  be 
a  “long-term”  r.v.,  i.e.,  it  is  constant  for  at  least  the  duration 
of  a  packet  slot. 

In  order  to  take  shadowing  into  account  in  the  analysis  of 
Section  III,  the  definition  of  the  constant  S  in  (5)  for  the  case 
of  channel-inversion  PC  must  be  modified  as  (5pc  =  7rr(l  — 
where  S  models  the  shadowing  between 
TXi  and  RXq  and  S'  models  the  shadowing  between  TX^ 
and  its  corresponding  RX.  In  the  absence  of  PC,  we  have 
i^npc  =  7rr(l  —  where  S'  now  denotes  the 

shadowing  between  TXq  and  RXq.  Therefore,  in  the  case  of 
no  PC,  Pi  in  (8)  is  conditioned  on  the  realization  of  S'  and 
the  expectation  of  the  former  with  respect  to  the  latter  must 
be  taken  in  order  to  obtain  the  unconditional  probability  of 
codeword/packet  error. 


The  sequence  {j/fe,  7fc,csi}^Jf  is  the  input  to  decoder  (6). 

A  consequence  of  channel  estimation  is  the  loss  of  infor¬ 
mation  rate  by  a  factor  MIM^,  due  to  the  transmission  of  the 
pilot  symbols  in  each  dwell.  This  factor  must  be  taken  into 
account  when  evaluating  the  transmission  capacity. 


B.  Correlation  of  the  interference 

Eq.  (8)  is  based  on  the  assumption  that  the  coded  bits  en¬ 
counter  independent  SIR  conditions  across  the  span  of  an  error 
event  in  the  decoder.  The  assumption  is  justified  if,  (a)  there  is 
a  sufficient  number  of  dwells^  such  that,  due  to  interleaving, 
these  bits  will  be  transmitted  on  different  dwells,  and  (b)  the 
SIRs  are  independent  across  dwells.  The  latter  assumption  is 
reasonable  if  the  number  of  frequencies  is  sufficiently  large. 
As  shown  in  [21],  the  temporal  correlation  of  the  interference 
in  a  fixed  Poisson  network  of  TXs  and  Rayleigh  fading  is  p/2, 
where  p  is  the  random  access  probability.  In  this  paper,  random 
access  is  achieved  via  random  FH,  i.e.,  p  can  be  substituted 
by  1/M. 


VI.  Numerical  results 

In  this  section,  we  present  numerical  results  for  a  network 
with  default  parameters  i?  =  1  m,  M  =  100,  Me  =  2  (QPSK), 
N  =  2  and  5  =  4.  Six  convolutional  codes  with  rate  i?c  =  1/2 
are  considered,  with  memory,  minimum  Hamming  distances 
and  distance  spectra  listed  in  Table  II  (see  [22]  for  more 
details).  In  all  plots,  the  upper  bound  to  the  pairwise  error 
probability  Pi  is  evaluated  by  (19). 

In  the  simulations,  we  generate  a  new  network  realization 
for  every  transmitted  packet.  The  network  area  is  selected 
such  that,  on  average,  there  are  60  interferers  per  subband.  At 
the  TX,  we  employ  a  block  interleaver  of  vertical  dimension 
equal  to  the  number  of  bits  per  dwell  MfTd  and  horizontal 


TABLE  II 

Optimum  rate  1/2  convolutional  codes 


Encoder 

Memory 

L 

Wl,  •  •  •  ,-!«L+4 

Enel 

1 

3 

1,1, 1,1,1 

Enc2 

2 

4 

2,0,5,0,13 

Enc3 

3 

5 

1,0,6,0,16 

Enc4 

3 

6 

1,3,5,11,25 

Enc5 

5 

8 

2,7,10,18,49 

Enc6 

7 

10 

1,6,13,20,64 

®This  is  also  related  to  the  length  of  the  codeword. 


Fig.  1 .  FEP  vs.  for  the  codes  in  Table  II.  The  solid  (dotted)  lines  depict 
(9)  with  error  events  up  to  length  L+4  (L)  taken  into  account.  The  theoretical 
result  is  slightly  optimistic  in  the  case  of  Enc6.  This  is  attributed  to  the  impact 
of  channel  correlation  within  the  packet,  the  effect  of  which  becomes  more 
apparent  as  the  span  of  the  error  events  becomes  larger  (QPSK,  L^,  =  800, 
D  =  40,  Ta  =  20,  N  =  2,  b  =  4,  perfect  CSI). 

dimension  equal  to  D,  such  that  D  consecutive  codeword 
bits  are  guaranteed  to  be  transmitted  on  different  dwells.  The 
coherence  time  of  the  fading  is  taken  to  be  equal  to  the 
duration  of  a  packet.  In  this  manner,  we  can  compare  the 
scheme  proposed  in  this  paper  with  the  slow-FH  (lA)  scheme 
advocated  in  [6], 

In  Fig.  1,  we  plot  the  FEP  vs.  N~^  (recall  that  A^eff  is  the 
expected  number  of  TXs  in  a  disc  of  radius  R  per  subband) 
for  the  codes  listed  in  Table  II  and  Lb  =  800.  The  upper 
bound  (9)  is  plotted  with  a  solid  line,  when  error  events  up  to 
length  L  +  4  are  taken  into  account,  and  with  a  dotted  line, 
when  only  the  Hamming  distance  error  events  are  taken  into 
account.  As  expected,  for  each  code,  the  two  curves  converge 
as  A^eff  — *■  0.  We  also  simulate  the  performance  of  the  codes 
by  dividing  the  packet  in  Z?  =  40  dwells  of  Td  =  20  QPSK 
symbols.  Fig.  1  demonstrates  that  the  dotted  curve  is  quite 
accurate  for  Pg  <  0.01. 

In  Fig.  2,  we  investigate  for  what  range  of  values  of  M  chan¬ 
nel  conditions  across  dwells  can  be  considered  independent. 
We  obtain  the  FEP  of  Enc4  via  simulation  for  different  systems 
which  have  the  same  A"eff  but  varying  M.  It  is  observed  that, 
for  M  >  20,  full  diversity  can  practically  be  attained,  as  the 
simulated  curves  are  very  close  to  the  theoretical  one  which 
is  derived  based  on  the  assumption  of  independent  channels. 
The  simulated  PEP  in  the  case  of  slow  PH  is  also  shown  for 
comparison.  Note  that,  given  Ngg,  the  performance  of  slow 
PH  does  not  depend  on  M;  its  diversity  order  is  always  one. 

In  Pig.  3,  we  plot  the  maximum  contention  density  A^, 
evaluated  by  (33),  vs.  L,  for  the  set  of  parameters  of  Pig.  1 
and  different  values  of  the  ratio  ejwL.  Note  that,  for  e/wL  = 
0.001,  increasing  the  diversity  order  of  the  code  from  3  to 
8,  results  in  a  tenfold  increase  of  A,;.  The  gain  comes  at  an 
increase  of  the  decoder  memory  from  1  to  5. 

Pig.  4  presents  the  results  of  Proposition  2  on  the  link 
ergodic  capacity  for  two  different  values  of  N .  As  expected, 
(25)  and  (26)  become  tight  as  A^eff  ^  oo  and  A"eff  ^  0, 


Ki  (dB) 


Fig.  2.  FEP  vs.  N~^  for  Enc4  and  M  =  4,  20, 100.  As  M  decreases  the 
diversity  order  decreases,  as  the  fading  and  interference  across  dwells  become 
more  correlated.  The  simulated  performance  for  slow  EH  is  also  shown  for 
comparison  (QPSK,  Lp  =  800,  D  =  40,  Td  =  20,  N  =  2,  b  =  4,  perfect 
CSI). 


Fig.  3.  Maximum  contention  density  vs.  L  for  different  values  of  ejwi^ 
(/?  =  1  m,  QPSK,  Lb  =  800,  N  =  2,  b  =  4,  perfect  CSI). 


respectively.  In  Pig.  5,  we  employ  the  channel  estimation 
scheme  of  Section  V  and  examine  the  effect  of  imperfect  CSI 
via  simulation,  when  Lb  =  500,  D  =  25  and  Td  =  20.  As  Tp 
is  increased  the  performance  loss  is  reduced,  at  the  expense  of 
information  rate,  e.g.,  at  Tp  =  10,  the  rate-loss  factor  is  2/3. 

Pinally,  in  Pig.  6,  the  simulated  PEP  of  Enc2  and  Enc4 
are  plotted  vs.  N^,  for  a  channel  with  lognormal  shadowing 
(as  =  6  dB)  and  shadowing-inversion  PC  or  no  PC.  Por  small 
A^eff,  PC  introduces  a  substantial  gain,  e.g.,  for  Enc4  at  N~g  = 
10  dB,  this  gain  is  an  order  of  magnitude.  The  reason  for  this 
is  that  the  detrimental  effect  of  shadowing  to  RXq  is  canceled 
and,  at  the  same  time,  the  interferer  process  is  sparse  enough 
such  that  the  code  protects  RXq  in  the  unlikely  event  of  a  large 
interfering  power  from  a  nearby  interferer  (this  occurs  when 
deep  shadowing  afflicts  the  channel  between  that  interferer  and 
its  respective  RX).  Theoretically,  we  can  see  this  with  the  help 
of  Jensen’s  inequality:  Given  that  all  shadowing  variables  are 
i.i.d.,  by  the  definitions  of  (5pc  and  (5npc  in  Section  V-D,  we 
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Fig.  4.  Link  ergodic  capacity  vs.  .  The  two  groups  of  curves  correspond 
to  TV  =  2, 10  (b  =  4). 


(dB) 


Fig.  5.  FEP  vs.  for  Enel  and  different  values  of  Tp.  For  Tp  =  5  and 
a  rate-loss  factor  of  20%  there  is  a  1  dB  penalty  compared  to  the  perfect-CSI 
case  (QPSK,  Lb  =  500,  D  =  25,  Ta  =  20,  N  =  2,  b  =  4). 

have  that  <  E[S'“']E[S'“]'  =  (5^p(,.  Hence, 

from  Proposition  1 ,  PC  should  perform  better  in  the  small  FEP 
regime,  which  is  verified  by  Fig.  6.  On  the  other  hand,  as 
the  network  becomes  very  dense,  the  diversity  in  the  received 
signal  is  lost  and  the  decoder  is  overwhelmed  by  an  increased 
interference  level  due  to  PC.  This  is  more  apparent  for  Enc2 
than  Enc4,  since  Enc2  has  a  smaller  Hamming  distance.  On 
a  final  note,  we  observe  that,  in  the  case  of  PC,  there  is  good 
agreement  between  the  theoretical  FEP  evaluated  by  (9)  and 
the  simulation  results. 

VII.  Concluding  remarks 

In  this  paper  we  considered  EH  during  packet  transmission 
and  coding,  as  a  physical-layer  scheme  for  random  wireless 
networks  with  uncoordinated  transmissions.  We  demonstrated 
via  analysis  and  simulation  that  the  transmission  capacity 
scales  as  where  e  is  the  constraint  placed  on  the  packet 
error  probability  and  L  is  the  code  diversity  order.  A  byproduct 
of  our  analysis  was  the  derivation  of  a  compact  expression 
for  the  pdf  of  the  SIR  in  a  Rayleigh  fading  and  a-stable 


Fig.  6.  FEP  vs.  for  Enc2  and  Enc4  with/without  PC  (QPSK,  Lb  =  500, 
D  =  25,  Td  =  20!  AT  =  2,  b  =  4,  perfect  CSI,  CTs  =  6). 

interference  channel,  when  the  RX  performs  MRC.  Upper  and 
lower  bounds  on  the  ergodic  capacity  of  this  channel  were  also 
derived. 

Employing  a  simple  channel  estimation  algorithm  based  on 
the  transmission  of  pilot  symbols  at  the  beginning  of  each 
dwell,  we  showed  that  the  performance  degradation  due  to 
imperfect  CSI  is  reasonable,  at  a  rate  loss  of  the  order  of 
20%.  The  effect  of  channel-inversion  PC  was  also  investigated 
for  a  channel  with  lognormal  shadowing,  and  PC  was  shown 
to  be  beneficial  if  A/M  is  sufficiently  small.  In  conclusion, 
we  believe  that,  given  the  gains  in  terms  of  network  capacity 
at  moderate  encoding/decoding  complexity,  even  for  a  small 
number  of  subbands,  this  scheme  merits  consideration  despite 
the  increased  overhead  compared  to  a  slow  EH  system. 
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